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Abstract 

In this paper we improve an earlier result by Bukhgeim and Uhlmann 1 , by showing that 
in dimension n > 3, the knowledge of the Cauchy data for the Schrodinger equation measured 
on possibly very small subsets of the boundary determines uniquely the potential. We follow 
the general strategy of but use a richer set of solutions to the Dirichlet problem. This 
implies a similar result for the problem of Electrical Impedance Tomography which consists in 
determining the conductivity of a body by making voltage and current measurements at the 
boundary. 
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1 Introduction 

The Electrical Impedance Tomography (EIT) inverse problem consists in determining the electrical 
conductivity of a body by making voltage and current measurements at the boundary of the body. 
Substantial progress has been made on this problem since Calderon's pioneer contribution [3], 
and is also known as Calderon's problem, in the case where the measurements are made on the 
whole boundary. This problem can be reduced to studying the Dirichlet-to-Neumann (DN) map 
associated to the Schrodinger equation. A key ingredient in several of the results is the construction 
of complex geometrical optics for the Schrodinger equation (see |14| for a survey). Approximate 



complex geometrical optics solutions for the Schrodinger equation concentrated near planes are 
constructed in [5] and concentrated near spheres in 

Much less is known if the DN map is only measured on part of the boundary. The only 
previous result that we are aware of, without assuming any a-priori condition on the potential 
besides being bounded, is in pQ. It is shown there that if we measure the DN map restricted 
to, roughly speaking, a slightly more than half of the boundary then one can determine uniquely 
the potential. The proof relies on a Carleman estimate with an exponential weight with a linear 
phase. The Carleman estimate can also be used to construct complex geometrical optics solutions 
for the Schrodinger equation. We are able in this paper to improve significantly on this result. We 
show that measuring the DN map on an arbitrary open subset of the boundary we can determine 
uniquely the potential. We do this by proving a more general Carleman estimate (Proposition 
3.2) with exponential non-linear weights. This Carleman estimate allows also to construct a much 
wider class of complex geometrical optics than previously known (section 4). We now state more 
precisely the main results. 

In the following, we let £1 CC R", be an open connected set with C°° boundary. For the 
main results, we will also assume that n > 3. If q € i°°(17), then we consider the operator 
—A + q : L 2 {Q) — ► L 2 (fl) with domain H 2 (tt) n Hy(il) as a bounded perturbation of minus the 
usual Dirichlet Laplacian. —A + q then has a discrete spectrum, and we assume 

is not an eigenvalue of - A + q : H 2 (il) n H^(ft) -> L 2 (n). (1.1) 

Under this assumption, we have a well-defined Dirichlet to Neumann map 

N q : JT*(0fi) d u u lm £ £T*(0fi), (1.2) 

where v denotes the exterior unit normal and u is the unique solution in 

H A (fl) := {u g H 1 ^); Au g L 2 (Q)} (1.3) 

of the problem 

(-A + q)u = in O, u\ m = v. (1.4) 

See pQ for more details, here we have slightly modified the choice of the Sobolev indices. 

Let Xq G R" \ ch (fi), where ch (O) denotes the convex hull of fi. Define the front and the back 
faces of dVl by 

F(x ) = {xe dVl; (x - x Q ) ■ v(x) < 0}, B(x ) = {x g dQ; (x - x ) ■ v(x) > 0}. (1.5) 
The main result of this work is the following: 

Theorem 1.1 With Q, xq, F(xo), B(xq) defined as above, let (71,(72 £ L°°(Q) be two potentials 
satisfying and assume that there exist open neighborhoods F, B C dfl of F(xq) and B(xq) U 

{x G dQ; (x — xq) ■ v = 0} respectively, such that 

N qi u = N q2 u in F, for all u e H?(dQ) n£'(B). (1.6) 

Then qi = q%. 

Notice that by Green's formula TV* = Nq. It follows that F and B can be permuted in (|1.6(l 
and we get the same conclusion. 

If B = d£l then we obtain the following result. 
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Theorem 1.2 With Q,, xq, F(xq), B(xq) defined as above, let 91,(72 € L°°(Q) be two potentials 
satisfying and assume that there exists a neighborhood F C dfl of F(xq), such that 

M qx U=N q2 v, in F, for all u G (1.7) 

Then q\ = qi- 

We have the following easy corollary, 

Corollary 1.3 With f2 as above, let x± G d£l be a point such that the tangent plane H of dfl at x\ 
satisfies <90 n H = {xi}. Assume in addition, that f2 is strongly starshaped with respect to x\. Let 
9i,Q2 G L°°(Q) and assume that there exists a neighborhood F C dft of x\, such that \1. 7| ) holds. 
Then q\ = qi- 

Here we say that Q is strongly star shaped with respect to x\ if every line through x\ which 
is not contained in the tangent plane H cuts the boundary dfl at precisely two distinct points, x\ 
and X2, and the intersection at xi is transversal. 

Theorem 1 1 . 1 1 has an immediate consequence for the Calderon problem. 

Let 7 e C 2 (f2) be a strictly positive function on fL Given a voltage potential / on the boundary, 
the equation for the potential in the interior, under the assumption of no sinks or sources of current 
in f2, is 

div(7Vit) = in f2, u \g n = /• 
The Dirichlet-to-Neumann map is defined in this case as follows: 

M,(J) = (fd v u) lm . 

It extends to a bounded map 

As a direct consequence of Theorem 1 . 1 we have 
Corollary 1.4 Let 7, G C 2 (fl), i = 1,2, be strictly positive. Assume that 71 = 72 on d£l and 
W 7l u = W 72 u in F, for all u G (c?Sl) n 

Then 71 =72- 

Here F and B are as in Theorem ll.il It is well known (see for instance ^3]) that one can relate 
7V 7 and J\f q in the case that q = with 7 > by the formula 

K(f) = (7-*)| fiI A(7-'/3 + ^(7- 1 M| 9n /- (1-8) 

The Kohn-Vogelius result 9 implies that 71 = 72 and d v 71 = 9„72 on F n B. Then using (|1.8|) 
and Theorem ll.ll we immediately get Corollary 1 1.41 

A brief outline of the paper is as follows. In section 2 we review the construction of weights 
that can be used in proving Carleman estimates. In section 3 we derive the Carleman estimate 
(Proposition 3.2) that we shall use in the construction of complex geometrical optics solutions 
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for the Schrodinger equation. In sections 4, 5 we use the Carleman estimate for solutions of the 
inhomogeneous Schrodinger equation vanishing on the boundary. This leads to show that, under 
the conditions of Theorems 1.1 and 1.2, the difference of the potentials is orthogonal in L 2 to a 
family of oscillating functions which are real-analytic. For simplicity we first prove Theorem 1.2. 
In section 6 we end the proof of Theorem 1.2 by choosing this family appropriately and using the 
wave front set version of Holmgren's uniqueness theorem. Finally in section 7 we prove the more 
general result Theorem 1.1. 

Acknowledgments. The first author was supported in part by NSF and at IAS by The von Neumann 
Fund, The Weyl Fund, The Oswald Veblen Fund and the Bell Companies Fellowship. The second 
author was partly supported by the MSRI in Berkeley and the last author was partly supported 
by NSF and a John Simon Guggenheim fellowship. 

2 Remarks about Carleman weights in the variable coeffi- 
cient case 

In this section we review the construction of weights that can be used in proving Carleman esti- 
mates. The discussion is a little more general than what will actually be needed, but much of the 
section can be skipped at the first reading and we will indicate where. 

Let C R", n > 2 be an open set, and let G(x) = (g % i{x)) a positive definite real symmetric 
n x n-matrix, depending smoothly on i e SI. Put 

p{x,0 = (G(x)£\0. (2.1) 

Let (f <G C°°(Q; R) with <p'(x) ^ everywhere, and consider 

p(x, £ + iift x {x)) = a{x, + ib(x, 0, (2.2) 

so that with the usual automatic summation convention: 

a(x, = g^ (x)tej - gV {x)tf/ Xi <p' x . , (2.3) 

b(x,0 = 2(G(x) l p'(x)\0 = 20"V*>- (2-4) 

Readers, who are not interested in routine calculations, may go directly to the conclusion of this 
section. 

A direct computation gives the Hamilton field H a — a'^ ■ d x — a' x ■ of a: 

H a = 2gV(x)£ j d Xi - d x unc^'k + ^(<7 ij 'KXA + 2<?V^XA' (2.5) 

and 

l -H a b = 2 5 %^>^> + 2 5 %^( 3 n^^ (2.6) 

= 2{^ x \Gi ® GO + 2(^ x \G<p' x ® Gip'J + 2(d x G\G£ ® V ' x ® 
— {d x G\G<p' x ® C ® + (dxG\G<pf x ® v' x ® <ff x ). 
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Here we use the straight forward scalar products between tensors of the same size (2 ohr 3) and 
consider that the first index in the 3 tensor d x G is the one corresponding to the differentiations 
d x . . We also notice that <p' x , £ are naturally cotangent vectors, while Gtp' x , G£ are tangent vectors. 
We want this Poisson bracket to be > or even e on the set a = b = 0, i.e. on the set given by 

(G|£ ® £ - *4 ® *4> = 0, (GK ®£>=0. (2.7) 

Observation 1. If ^ is a distance function in the sense that (G\ip' x <E> f' x ) = 1, then if we differentiate 
in the direction G(p' x , we get 

= (G<p' x ■ d x )(G\<p x <g> <p' x ) = {d x G\G v ' x ®v' x ® V' x ) + 2(V XX \G^ X <8> G^>. 
From this we see that two terms in the final expression in l|2.(j[l cancel and we get 

~H a b = 2(^\G£ ® GO + 2(d x G\G^ ® <f/ x ® - (d x G\G<p' x ® £ ® £>■ ( 2 -8) 



Observation 2. If we replace y(x) by V'(x) = f(tf(x)), then 

< = fMxW x 

= />(xm®^+/v*m*. 

If £ satisfies (I2.7|) . then it is natural to replace £ by r\ = /'(y>)£, in order to preserve this condition 
(for the new symbol) and we see that all terms in the final member of l|2. 6fl . when restricted to 
a = b = 0, become multiplied by /'(</?) 3 except the second one which becomes replaced by 

/'(^) 3 2(^JG^ ® G^) + 2/Xz))/' (^)) 2 (G|^ ® ^) 2 - 

(For the first term in we also use that ® ^|G£ ® G£) = |G£) 2 = 0.) Thus we get after 
the two substitutions tp\-* ij) = f(ip(x)), £ 1— > 77 = f'(cp(x))£: 

\H a b{x,ri) = 2f"{ V {x))f{ V {x)Y\WX+ (2-9) 
/'M 3 (2(^'JG£ ® GO + 2< ¥ £ X |<V B ® G^> + 2(a c G|G£ ® ^ ® 
—{d x G\G(p' x ® £ ® + ( 9 xG\G(p' x ® ^ ® ip' x )), 

with 7] = /'(</?)£, £ satisfying (|2.7I) . so that 77 satisfies the same condition (with replaced by -0): 

(G^ (8 t? - < « O = (G|< ® //) = 0. (2.10) 

Moreover ||^||g = {G\(f' x ® <^4) by definition. 

Conclusion. To get i? a 6 > whenever (|2.7|) is satisfied, it suffices to start with a function <p 
with non-vanishing gradient, and then replace ip by f(ip) with /' > and /"//' sufficiently large. 
This kind of convexification ideas are very old and used recently in a related context by Lebeau- 
Robbiano JEj, Burq 0. For later use, we needed to spell out the calculations quite explicitly. 
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3 Carleman estimate 

We use from now on semiclassical notation (see for instance 0]). 

Let Po = —h 2 A = ^(/iDr.) 2 , with D x . = id Xj . Let ip, fi be as in the beginning of Section|21 
Then 

n 

e v ' h o P o e-*' h = ^2(hD Xj + id Xj pf = A + iB, (3.1) 

where A, B are the formally selfadjoint operators: 

A = {hDf - {p' x f, B = V + hD Xj o X .^) (3.2) 

having the Weyl symbols (for the semi-classical quantization) 

a = e-{v' x ?, b = 2<pl,>Z. (3.3) 

We assume that tp has non vanishing gradient and is a limiting Carleman weight in the sense 
that 

{a, b}(x, £) = 0, when a(x, £) = b(x, f) = 0. (3.4) 
Here {a, b} — a'^ ■ b' x — a' x ■ b'^ is the Poisson bracket (as in (|2.tj[l 'l: 

{a,6} = 4<^ x (x)|£®£ + <^®^>. (3.5) 

On the rr-dependent hypersurface in £-space, given by b(x,^) = 0, we know that the quadratic 
polynomial {a,b}(x,£,) vanishes when £ 2 = (p' x ) 2 - It follows that 

{a, b}(x, = c(x)(e - (^) 2 ), for b(x, - 0, (3.6) 

where c(x) € C°°(f2;R). Then consider 

{0,6X1,0 -cCa:)^ 2 -^) 2 ), 

which is a quadratic polynomial in £, vanishing when </? x (x) • £ = It follows that this is of the 
form £(x, £)6(x, £) where £(x, £) is affine in £ with smooth coefficients, and we end up with 

{a, b} = c(x)a(x, + £(1, 06(z, £)• (3-7) 

But {a, 6} contains no linear terms in £, so we know that £(x, £) is linear in £. 
The commutator [A, S] can be computed directly: and we get 

= \{Ytt hD *i ° ohD X] )hD Xk +hD Xk (hD X] o^. Xk +pl ]Xh ohD X] )\ 

+4<V^ a ,¥4(as)®¥4(s))). 
The Weyl symbol of [A, B] as a semi-classical operator is 

-{a,b} + h 3 p (x), 
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Combining this with (|3.7[) . we get with a new pq: 

i[A, B] = h(^(c(x) oA + Aoc) + - (LB + BL) + h 2 p {x)), (3.8) 

where L denotes the Weyl quantization of I. 

We next derive the Carleman estimate for u £ Cq°(0), fl CC £1: Start from P$u = v and let 

u = e v / h u, v = efl h v, so that 

{A + iB)u = v. (3.9) 
Using the formal selfadjointness of A, B, we get 

||v|| 2 = ((A - iB)(A + iB)u\u) = \\Au\\ 2 + \\Bu\\ 2 + (i[A, B]u\u). (3.10) 

Using EHJ), we get for u € Cg°(Q): 

\\v\\ 2 > ^ II^^TII 2 <^C^) Cll^^ll ll^ll ll^^ll ) ^C^ 3 )II^H 2 (3.H) 

> ^\\Auf + ^\\Bu\\ 2 -0(h 2 )(\\u\\ 2 + \\Luf). 

> i(||^|| 2 + ||^|| 2 )-0(^)P|| 2 , 

where in the last step we used the apriori estimate 

||W£|| 2 <0(1)(||^|| 2 + H| 2 ), 

which follows from the classical ellipticity of A. 

Now we could try to use that B is associated to a non-vanishing gradient field (and hence 
without any closed or even trapped trajectories in f2), to obtain the Poincare estimate: 

h\\u\\<0{l)\\Bu\\. (3.12) 

We see that <|3.12ll is not quite good enough to absorb the last term in (|3.11(l . In order to 
remedy for this, we make a slight modification of ip by introducing 

fe = f o <p, with / = f e (3.13) 

to be chosen below, and write a e +ib e for the conjugated symbol. We saw in Section|21and especially 
in (|2.9|) that the Poisson bracket {a e ,b e }, becomes with <p equal to the original weight: 

K, b t }(x, f\<p)rj) = f'(v) 3 ({a, b}(x, n) + y^ll^H 4 )' when V) = V) = 0- (3.14) 

The substitution £ — > f'(<p)r) is motivated be the fact that if a(x, rj) = b(x, rj) — 0, then a e (x, f{<p)rj) - 
b e (x,f'(<p)i}) = 0. Now let 

MX) =A + eA 2 /2, (3.15) 



with < £ < 1, so that 

4£ -±e + 0(e 2 ). 



f'(<P) 1 + ^ 
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In view of [ETHfl . (|3~3|) . we get 

{a e ,b e }(x,0 = 4/ e »(/^)) 2 || ¥ /|| 4 « 4 e ||^|| 4 , (3.16) 
when a e (a;,£) = 6 c (x,£) = 0, so instead of (|3.7|) . we get 

{a £ ,6 £ } = 4/ £ "(^)(/ e '(^)) 2 ||^|| 4 + c e (x)a e (x,C)+4(x,0& £ (^,0, (3-17) 

with £ € (o;,£) linear in £. 

Instead of (|3.11|) . we get with u = e v ^ h u, v — e v '^ h v when P u = v: 

\\v\\ 2 > h{Ae + 0(e 2 )) J yj*\u(x)\ 2 dx + ±\\A e u\\ 2 + ±\\B t u\\ 2 C(^ 2 )ll"l| 2 , (3.18) 

while the analogue of l|3.12[l remains uniformly valid when e is small: 

h\\u\\ < 0(l)\\B e u\\, (3.19) 

even though we will not use this estimate. 
Choose /iC e< 1, so that (|3.18|l gives 

\\v\\ 2 >eh\\u\\ 2 + l\\A e u\\ 2 + ±\\B e u\\ 2 . (3.20) 

We want to transform this into an estimate for u 1 v. From the special form of A e , we see that 

\\hDu\\ 2 < (A e u\u) + 0(l)\\u\\ 2 , 

leading to 

\\hDu\\ 2 < ^\\A e u\\ 2 + 0{l)\\u\\ 2 . 
Combining this with (|3.2U|) . we get 

\\vf > ^(Pll 2 + \\hDuf) + (1 - 0{A))\\ASf + l\\B e u\\ 2 . (3.21) 

Write ip e = Lp + eg, where g = g t is 0(1) with all its derivatives. We have 

u = e eB ' h u, v = e t9,h v, 

so 

e 



hDu = e t9,h {hDu + -g'u) = e^ 9 (hDu + 0(e)u), 
i 

and 

||w|| 2 + \\hDu\\ 2 > \\e C9 ' h u\\ 2 + \\e t9 / h hDu\\ 2 - Ce\\e t9 ' h u\\ \\e t9/h hDu\\ - Ce 2 \\e £g/h u\\ 2 
> (1 - Ce){\\e t9 ' h u\\ 2 + \\e e9 ' h hDu\\ 2 ), 

so from (|3.21|l we obtain after increasing Cq by a factor (1 + 0(e)): 

\\e t9/h v\\ 2 > ^{\\e t9/h u\\ 2 + \\e e9 / h hDu\\ 2 ). (3.22) 
L-0 
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If we take e = Ch with C» 1 but fixed, then eg/h is uniformly bounded in S7 and we get the 
Carleman estimate 

h 2 (\\u\\ 2 + \\hDu\\ 2 ) < dWvf. (3.23) 
This clearly extends to solutions of the equation 

(-h 2 A + h 2 q)u = v, (3.24) 

if q € L°° is fixed, since we can start by applying (|3.23|) with v replaced by v — h 2 qu. Summing 
up the discussion so far, we have 

Proposition 3.1 Let Pq, Q, tp be as in the beginning of this section and assume that (p is a limiting 
Carleman weight in the sense that \S.4\j holds. Let Q, CC O be open and let q € L°°(S1). Then if 
u e Cq?(Q), we have 

h(\\e v / h u\\ + \\hDe*/ h u\\) < C\\e^ h (-h 2 A + h 2 q)u\\, (3.25) 

where C depends on SI, and h > is small enough so that C/i||q||£=o(Q) < 1/2. 

We next establish a Carleman estimate when Pqu — v, u G (7°°(fi), W| ao = and SI CC ft is 

a domain with C°° boundary. As before, we let u = e v ^ h u, v = e v ^ h v, with ip = (p e: < e^C 1. 
With A = A e , B = B e , we have 

{A + iB)u = v, (3.26) 

and 

||u|| 2 = {{A + iB)u\{A + iB)u) (3.27) 
= ||Au|| 2 + \\Bu\\ 2 + i((Bu\Au) - (Au\Bu)), 

Using that B is a first order differential operator and that 

«lan = °> 

we see that 



(Au\Bu) = (BAu\u). (3.28) 

Similarly, we have 

(Bu\( V ' x ) 2 u) = ((p' x ) 2 Bu\u). (3.29) 
Finally, we use Green's formula, with v denoting the exterior unit normal, to transform 

{Bu\ - h 2 Au) n = -h 2 (Bu\d v u) dn + (-h 2 ABu\u) n , 

where we also used that u\ dn = 0. 
On <9S7, we have 

B = 2(<p' x -v)~d v +B', 
i 

where B' acts along the boundary, so using again the Dirichlet condition, we get 

(Bu\d v u) d n = —{{<p' x ■ v)d v u\d y u) d Q,. 
i 
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Putting together the calculations and using l|3.2|) for A, we get 

||v|| 2 = \\Au\\ 2 + \\Bu\\ 2 + i([A, B]u\u) - 2h 3 {{^' x ■ u)d u u\d v u) aQ . (3.30) 

Let 

dn ± ={xe Oil; ±(p' x ■ v > 0}. (3.31) 
Notice that dfl± are independent of e. We rewrite (13.301) as 

-2h 3 ((^' x ^)d^d,u) dn _+i{[A,B}u\u) + \\A^ + \\Bu\\ 2 = \\vf+2h 3 (( i p' x -u)d l/ u\d l/ u) an+ . (3.32) 

This is analogous to 13.10fl and the extra boundary terms can be added in the discussion leading 
from l|3.18|l to l|3.21|l and we get instead of (I3.21|) : 

-2h 3 {{p' x ■ v)d v u\d v -u) d n_ + i|(N| 2 + \\hDu\\ 2 ) + (3.33) 

Cz - 0(eh))\\A £ u\\ 2 + h\\ B M 2 

< \\d\\ 2 + 2h 3 {{tp' x ■ v)d v u\d v u) aa+ , 

with ip — ip e , provided e~^> h. Fixing e = Ch for C > 1, we get with ip — <p e= o for some Co > 0: 

' v)d v u\d v u) dn _ + g(N| 2 + \\hDuf) (3.34) 

< \\v\\ 2 + C h 3 ((<p x ■ v)d u u\d„u) an+ . 

Here we recall that —h 2 Au = v, u = e v ^ h u, v = e v ^ h v, tp = (p e= o, u\ gsi = 0. 

If q G L°°, we get for h 2 (—A + q)u = v, it| an = 0, by applying l|3.34|l with v replaced by 
v — h 2 qu: 

• v)d v u\d v u) aa _ + f o (\\u\\ 2 + \\hDu\\ 2 ) (3.35) 

< ||v|| 2 + C a h 3 ((tp' x • v)d v u\d v u)an + . 

Here u, v are defined as before. 
Summing up, we have 

Proposition 3.2 Let f2, ip be as in Provosition \3.1\ Let fl CC be an open set with C°° boundary 
and let q £ L°°(fl). Let v denote the exterior unit normal to dfl and define dfl± as in \S.S1\) . 
Then there exists a constant Cq > 0, such that for every u S C°°(fl) with u\ g ^ = 0, we have for 
0<h<l: 

• v)e^ h d v u\e^ h d v u) dn _ + ^.(\\ e ^ h u\\ 2 + He^V^ 2 ) (3.36) 
< \\ev/ h {-h 2 A + h 2 q)u\\ 2 + C h 3 {{p' x ■ v)e*/ h d„u\e*/ h d u u) dn+ , 



Remark. If ip is a limiting Carleman weight, then so is —ip. With u = e~ v / h u, v = e~ v / h v, we still 
have (|3.35|) . provided we permute <9f2_ and dfl+ and change the signs in front of the boundary 
terms, so that they remain positive. 
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4 Construction of complex geometrical optics solutions 

Let 7J s (R ra ) denote the semi-classical Sobolev space of order s, equipped with the norm || (hD) s u\\ . 
We define 7J s (ft), H§(Q) in the usual way, when li CC R ra has smooth boundary. 1(3.23(1 can be 
written 

h\\u\\ m < C\\e v / h P e- v/h ul u £ Cg°(fi), (4.1) 

when P Q = -h 2 A. Here we let ft C ft be as in Sectional Recall that P a ^ = e v ^ h P Q e~ v/h has the 
semiclassical Weyl symbol £ 2 — tp' x 2 + 2iip' x ■ £ = a + ib, which is elliptic in the region |£| > 2|<p'(x)|. 
It is therefore clear that (|4.1|) can be extended to: 

h\\u\\ H -s +1 < C Si a\\e^ h P e-^' h u\\ H -s, u £ C °°(ft), (4.2) 

for every fixed s £ R. With q £ L°° (ft) , we put 

P = -h 2 (A -q), P v = e ^ h Pe~^ h = P Q . V + h 2 q. 

If < s < 1, we have 

IMIff- < IIH < NMMI < ||g|U-||u||ir--H-i 1 
and for h > small enough, we get from 1(4.2(1 : 

^a.nlle^Pe-^ullff-.. (4.3) 
The Hahn-Banach theorem now implies in the usual way: 

Proposition 4.1 Let < s < 1. TTien /or ft > small enough, for every v £ iP 7 (ft), tftere 
exists u £ iP(ft) such that 

e- lp/h Pe lp/h u = v, h\\u\\ a ><C\\v\\ B .-i. (4.4) 

This result remains valid, when o: is complex valued. In that case we replace P in 14.3(1 by 
P = -h 2 A + q. 

We next construct certain WKB-solutions to the homogeneous equation. Recall that a, b are 
in involution on the joint zero set J : a = 6 = in view of 1(3.7(1 . At the points of J we also see 
that the Hamilton fields 

H a = 2(£ ■ d x + (^VJ%», # b = 2(^4 • 5, - (<^|%» (4.5) 

are linearly independent and even have linearly independent x-space projections. We conclude 
that J is an involutive manifold such that each bicharacteristic leaf (of dimension 2) has a base 
space projection which is also a nice submanifold of dimension 2. It follows that we have plenty of 
smooth local solutions to the Hamilton- Jacobi problem 

a(x, ip'{x)) = b(x, ijj{x)) = 0. (4.6) 

Indeed, if (xo,£o) e J-> anc ^ we let P C ft be a submanifold of codimension 2 passing through xo 
transversally to the projection of the bicharacteristic leaf through (xo, £o), then we have a unique 
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local solution of l|4.6|) . with tp\ H = tp, if tp is a smooth real- valued function on iJ such that tp'(xo) 
is equal to the projection of (xo,£o) hi T* Q (H). 

Since we need some explicit control of the size of the domain of definition of tp, we now give a 
more down-to-earth construction. Ij4.6|l can be written more explicitly as 

tp'(x) 2 - ip'(x) 2 = 0, <p'(x) ■ tp'(x) = 0. (4.7) 

First restrict the attention to the hypersurface G = (y5 _1 (Co) for some fixed constant Co, and let g 
denote the restriction of tp to G. Then we get the necessary condition that 

g'(xf = V '{x)\ (4.8) 

where g'(x) 2 is the square of the norm of the differential for the metric dual to eo, the induced 
Euclidean metric. Now i|4.8|) is a standard eikonal equation on G and we can find solutions of 
the form g(x) = dist (x,T), where T is either a point or a hypersurface in G and dist denotes the 
distance on G with respect to the metric <p' (x) 2 eo(dx) . Of course, we will have to be careful, since 
such distance functions in general will develop singularities, and in the following we restrict G if 
necessary, so that the function g is smooth. With g solving 1)4. 8|) . we define tp to be the extension 
of g which is constant along the integral curves of the field f'{x) ■ d x : 

<f/(x)-d x 1>(x)=0, ip\ G = 9- (4.9) 

Then the second equation in (|4.7|l holds by construction, and the first equation is fulfilled at the 
points of G. In order to verify that equation also away from G, we consider, 

y'{x) ■ W 2 - v' 2 ) = m"v'W) - <¥>»'»• (4-10) 
Taking the gradient of the second equation in (|4.7|l . we get ip"ip' + ip"ip' = 0, and hence 

p'(x)-dM /2 -^ 2 ) - -2({</W) + WW)) = -l{a,b}(*rf) (4- 11 ) 

= -lc{x){^-ip'*)-l[x^'W -i,' 

= -l c (x)W 2 -<p> 2 ). (4.12) 

Thus 

fo/(s) • d x + - </) = 0, (^' 2 - ^' 2 )| G = 0, (4.13) 

and we conclude that tp' 2 — f'(x) 2 = 0. 

Summing up the discussion so far, we have seen that if tp is a limiting Carleman weight, and 
the open set U, is a union of integral segments of (p'(x) ■ d x all crossing the smooth hypersurface 
G C </3 _1 (Co), then if g is smooth solution to the eikonal equation \4-$ on G and we define tp to 
be the solution of {4-^ , we get a solution of l4-6\ ). 

(|4.(i|) implies that 

p(x, i(p'(x) + tp'(x)) = 0, (4.14) 
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which is the eikonal equation for the construction of WKB-solutions of the form u(x; h) — a(x; h)ek(~ 
of Pou 0. If we try a smooth and independent of h, we get 

= {((hD + 4,' x f - + i(<f/(x)(hD + ^) + (hD + tf)<ff)))a 
= (hL - h 2 A)a, 

where L is the transport operator given by 

L = iP'D + DiP' +i(<p'D + D<p'). (4.16) 

Along the projection of each bicharacteristic leaf this is an elliptic operator of Cauchy-Riemann 
type and if we assume that the leaves are open and simply connected, then (see 5 ) there exists a 
non-vanishing smooth function a S C°° such that 

La = 0. (4.17) 

Recall that q £ L°°(f2). Assume that a in 14. 1711 is well-defined in a neighborhood of f2. Then 
from H4.15fl . we see that with P = Pq + h 2 q: 

p e i(-<p+ii>) a = e -v/hh 2 d, (4.18) 

with d = 0(1) in L°° and hence in L 2 . Now apply Proposition 14. II with <p replaced by — <p, to find 
r € H 1 ^) with h\\r\\ H i < Ch 2 , such that 

eV /hp e - V /h e i^/h r _ _ h 2^ 

i.e. 

P( e K(-^+#)( a + r )) =o. (4.19) 

5 More use of the Carleman estimate 

In Section|31we derived a Carleman estimate for e lf ^ h u when h 2 (— A + g)tt = u when y> is a smooth 
limiting Carleman weight with non- vanishing gradient. In order to stick close to the paper pQ, we 
write the corresponding estimate for e~ tp ' h u, when (—A + q)u = v, U\ dQ = 0: 

• u)e-^ h d v u\e-^ h d u u) dn+ + g(||e-^u|| 2 + \\e^/ h hVu\\ 2 ) (5.1) 
< h A \\e-^/ h v\\ 2 - C h 3 ((p' x ■ v)e-^/ h d v u\e-f/ h d l ,u) d n_: 

where v is the exterior unit normal and Q± = {x € dQ] ±is ■ ip' > 0}. 
Let qi,q 2 G L°°(ri) be two potentials. Let 

u 2 = er^+^\a 2 + r 2 (x- h)), with (A - q 2 )u 2 = 0, ||r 2 || H i = 0{h). (5.2) 

Here -02 is chosen as in Section 0] so that (ip') 2 = (V4) 2 = f' • ip 2 = an d so that the integral 
leaves of the commuting vector fields tp' • d x ,tp' • d x are simply connected in Q. a 2 is smooth in a 
neighborhood of Q and everywhere non- vanishing. 
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Let N q be the Dirichlct to Neumann map for the potential q and let 

<9fi_. eo = {x S dil; v(x) ■ (f' x (x) < e }, 
dCl +t e = {x e dQ,; v{x) ■ <p' x (x) > e }, 

for some fixed eo > 0, so that dfl + ^ C <9f2_ C c?Sl_ :£o . Here v(x) denotes the unit outer 

normal to dft. 
Assume 

=Nq 2 (f), in £>n_, eo , for all / e fT*(0fi). (5.3) 

Let ux S J H" 1 (fi) solve 

(A-« 1 )« 1 =0, « 1 | 8n =« 2 | 8n . (5.4) 
Then by the assumption H5.3() . we have 

d v ui = d v u 2 in 90_ i£o . (5.5) 

Put u = ui — «2, q = 52 — gi, so that 

supp (di/U| an ) C 9fi+, eo , (5.6) 

and 

(A - g x )u = (A - qi)u 2 = qu 2 , u\ gQ = 0. (5.7) 
For u G 7J 1 (fi) with Aw £ L 2 (il), we get using (|5.6fl . l]5.7fl and Green's formula: 

/ qu 2 vdx = / (A — qi)uvdx = I u(A — q\)vdx + f (d v u)v S (dx) . (5-8) 
Jn Jn Jn Jao +j£0 

As in Section 21 we can construct 

v = e-^ +i ^(ai+n), (5.9) 

with ipi satisfying ip' ■ tpi — 0, {p 1 ) 2 — {ip'i) 2 , with a±(x) non-vanishing and smooth, and with 
lki||fl-i(n) = O(h), so that 

(A-q 1 )v = 0. (5.10) 

Then (|5.8|l becomes 

qe^ 1+ ^\a 2 + r 2 )(ai + n)dx = [ {d^e^-^ (ai + n) S(dx). (5.11) 
We shall work with ipi,ip 2 ,ip slightly /i-dependent in such a way that 

h*l>i+rh)^f, h ->0. (5.12) 
a 

Recall that 

IN| H i=0(fc). (5.13) 
Then using that g S we see that the left hand side of l|5.11|l converges to 

a 2 aiq(x)e lf{x) dx. (5.14) 
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For the right hand side of i|5.11fl . we have, using (|5.1|) . for (A — qi) and l|5.7|) : 



| f {d v u)e~^- l ^\ ai + n )S(dx)\ 2 



(5.15) 



< 



l|oi + ri||| n+> . / {e~*l h \d v u\) 2 S{dx) 



< 



\\ai+n\\l a ± [ {<p>-v){e-v' h \d v u\fS(dx) 



£ o Van 



< 




Here d v u = on c?fl_, and using also (|5.2ll . we get 



| / (d u u)e-^-^(ET+^)S(dx)\ 2 < ^||oi + n||L + ||g(a a + r 2 )| 




(5.16) 



Here ||<?(a2 + r 2 )\\ 2 = 0(1), by (|5.13|) . Since n = 0(h) in the semiclassical JT 1 -norm, we have 
r\ = 0(1) in the standard (h = 1) ii^-norm. Hence 



Consequently, the right hand side of l|5.11(l tends to 0, when h — ► 0, and letting h — > there, we 
get 



./n 

for all / that can be attained as limits in i|5.12[l . 

Finally, we remark that if ip is real-analytic, then in the above constructions, we may arrange 
so that tpj and a,j have the same property. 

6 End of the proof of Theorem 1.2 

From now on, we assume that the dimension n is > 3. We choose tp(x) = hi\x — xq\ for xq 
varying in a small open set separated from by some fixed affine hyperplane H . Notice that ip is 
a limiting Carleman weight in the sense of (|3.4Jl . We need a sufficiently rich family of functions / 
in (|5.18|) and recall that these functions are the ones that appear in l|5.12[l with ipj analytic near 
f2 and satisfying (ipj) 2 — (<p') 2 , ip'j • <p' = 0- Changing the sign of ipi we can also view / as a limit 
h(ipi — 1P2) for suitable such /i-dependent functions ipj. More precisely, we can take an analytic 
family ip(x, a) depending on the additional parameters a = (oti, atk), with ip(-, a) satisfying 



0(1) in L 2 . 



(5.17) 




(5.18) 



(V4) 2 = (^) 2 , <Vx = o, 



(6.1) 



and then take 




(6.2) 
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We first discuss the choice of ip. Since ip' x is radial, with respect to xq, the second condition in 
means that ip(x) is positively homogeneous of degree with respect to x — xq. A necessary 
and sufficient condition for ip (at least if we work in some cone with vertex at xq) is then that 

= W x )\ (6-3) 

on a suitable open subset xq + tqW of xq + roS 71 ^ 1 , for some fixed ro > 0. The necessity is obvious 
and the sufficiency follows easily by extending ip to be a positively homogeneous function of degree 
in the variables x — xq. 

Here is an explicit choice of a suitable open set in i|6.3[l : Let ro > be large enough so that 
O C B(x ,r ). Let x + r W C dB(x Q ,r ) be defined by 

x + r W = dB{x ,r )nH+, (6.4) 

where H + is the open half-space delimited by the affine hyper-plane H, for which xq (jt H + (so 
that Vt C H+). Then Q, is contained in the open cone xq + R+VT, so if we choose ip on xq + r$W 
as in (|6.3|) and extend by homogeneity, we know that ip will be smooth near Q. 

Let ?/o € dB(0, 1) \ W be such that the antipodal point — yo also is outside W and define 

ip(x,y) = d S n-i(x,y). (6.5) 

Then ip £ C°°(W x neigh (yo)) and the function ip((x — xq)/\x — xo\, y) € C°°(Cl x neigh (yo)) will 
satisfy (|6.1|l . Since the domain of definition does not contain antipodal points, we remark that 

is of rank n - 2 and ^«^) = W a ) x M^, v ) = «)• (6-6) 

This follows from basic properties of the geodesic flow (and remains true more generally for 
ip{x, y) = d(x, y) on a Riemannian manifold as long as x, y are not conjugate points.) 
For x e W C S n -\ (y, v) E TS^ 1 , y € neigh (y ), we put 

f(x;y,v)=ip' y (x,y)-u. (6.7) 

Then 

f x {x;y,v)=^ y {x,y){v). (6.8) 

In view of (|6.fi|l . we see that this vanishes precisely when v || ip' y (x, y), i.e. when v is parallel to the 
(arrival) direction of the minimal geodesic from x to y. Restricting v to non-vanishing directions 
which are close to be parallel to the plane H, we can assure that 

%(x;y,v)?0. (6.9) 

Lemma 6.1 f" , has maximal rank n — 1. 

J xAv.u) 



Proof: We already know that f'J v = ip 1 ' is of rank n — 2 and that the image of this matrix is 
equal to (tp' x ) ± . Consequently, we consider 

9(y) = ^' x (x,Vo) ■ ipx, y (x,y)(v) = {ip'J. ty (x,y)\ip' x (x,y a ) O v) 
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as a function of y £ neigh (yo). The function vanishes for y — yo and can also be written 

yW x {x,y ) - ^(x, y)) ®v) = (^ ty (x,y)(u)\^ y (x,y)(y - y)) + O((y - yf). 

From this expression, we see that the y-differential is non- vanishing and hence the range of /" / v s 
contains vectors that are not orthogonal to ^' x (x,y). □ 

Now consider 

*(x;y,x) = i>{^—^,y) £ C°° (H x neigh (y , S^ 1 ) x neigh (x 0) R n )). (6.10) 
\x — x\ 

^> is analytic, real and satisfies l|6.1[l with (p(x) = &(x, x) = In |x — x|. We can take a = y and (|6.2|) 
becomes 

f{x) = f{x;6) = {*' v \v), = {y,x,u), (6.11) 

with (y, v) £ TS 11 ^ 1 . Lemma 16.11 shows that f' x ^ v \ has rank n — 1 and indeed the image of 
this matrix is the tangent space of dB(x, \x — x\) at x. Since f x is a non-vanishing element of 
T x (dB(x, \x — x|)), we can vary x infinitesimally to see that f xx (Ji) ^ T x (dB(x, \x — x|)) for a 
suitable Jl £ R". It is then clear that 

f'x.e = f",{ y ,x,v) nas maximal rank n, (6.12) 

and hence that the map 

neigh (fi) 3 x i-» /£(x, 0) G R 3n ~ 2 (6.13) 

has injective differential. 

Lemma 6.2 TTie map ffi.l!])) is injective. 

Proof: Let Xi, X2 £ neigh (fl) be two points with 

fe(?i,9) = f'g(x 2 ,9), (6.14) 
for some # = (y, z/). Taking the ^-component of this relation, we get 

x ■ — x 

i>y(Vl,V) =i>'y{X2,y), Xj = y. (6.15) 

\Xj - x\ 

This means that x\, X2, y belong to the same geodesic 7 and this geodesic is minimal (i.e. 
distance minimizing) on some segment that contains these three points in its interior. If x\ 7^ X2, 
we may assume that d(x2, y) < d(xi, y). For y £ neigh (yo, S" 1-1 ), we have 

d(x 1 ,x 2 ) + d(x 2 ,y) - d{x x ,y) =: g(y), g(y) ~ d(y,j) 2 . 

It follows that 

/(x 2 ; y, x, v) - f(xi,y, x, v) = g'{y) ■ u, 

and using that v is not parallel to 7 at yo, we see that this function has a non- vanishing y-gradient 
at y = yo, in contradiction with H6.14[l . Thus, x\ = x 2 , or in other words, Xi and X2 belong to the 
same half-ray through x. 
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Taking the x-component of (|6.14[) . we get 

\x—X\ \x—X2 

These quantities are clearly non- vanishing and if x\ 7^ x 2 , they differ by a factor 7^ 1, since 
Vg( ) is homogeneous of degree — 1 in x — x. Thus x\ — x 2 . □ 

Now apply CT| with /(a?) = /(x,0): 

e if ^a 2 a iq (x)dx = 0, (6.16) 

where 02, ai are analytic non-vanishing functions of x, y, x in a neighborhood of f2 x {yo} x {^o}- 
Since f(x, 9) = f(x; y,x, v) depends linearly on v, we can replace v by \v and get 

z lXf{x > e) a 2 a iq (x)dx = 0, A > 1. (6.17) 

Now represent 9 by some analytic real coordinates 9i,9 2l near some fixed given point 9q = 

(yo, xq, vq). If x, z € fi, w € neigh (6q), we consider the function 

9^-f(z,9) + f(x,9)+ l -(9-w) 2 . (6.18) 

For x = z, we have the unique non-degenerate critical point 9 = w, while for i/z there is no real 
critical point in view of Lemma 16.21 For x ~ z we have a unique complex critical point which is 
close to W, and we introduce the corresponding critical value 

f/>(z,x,w) = v. C . e {-f{z,6)+f{x,6) + ~{6-w) 2 ). (6.19) 

From (|6.13|) and standard estimates on critical values in connection with the complex stationary 
phase method (^DE3)> we deduce that 

lmi/j(z, x, w) ~ (z — x) 2 , z,x£fl, z w x, (6.20) 

Moreover, when x = z, we have 

t/> z (z,z,w) = -f z (z,w), ip' x (z,z,w) = f z (z,w), ip(z,z,w) = 0. (6.21) 

We now multiply (|6.17|l by x(9 ~ w } e lX %( e - w ) 2 - lX f( z , e ) j and integrate with respect to 9, to get 

e lX ^ z ' x ' w) a(z, x, w; X) X (z - x)q(x)dx = 0{e~^). (6.22) 



Here X denote (different) standard cutoffs to a neighborhood of 0, and a is an elliptic classical 
analytic symbol of order 0. 

Now restrict w to an rt-dimcnsional manifold S which passes through 6*o, and write (z, —f' z (z, 9)) = 
(a x , a^) — a. Then wc rewrite ()6.22(l as 



J e ^(^) a ^ x] h) x (a x -x)q{x)dx =0(e~^), 



(6.23) 
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implying that 

(z,-f' z (z,9 ))^WF a (q), (6.24) 
since we can apply the standard FBI-approach (^3]). Notice that l|6.20[l . (|6.21() give: 

tp(a, x) — (a x — x) ■ + 0((a x — x) 2 ), Imip(a, x) ~ {a x — x) 2 , (6.25) 

and we can choose £ so that the map neigh (zq) xE 9 (z, ff) i— > (z, —f' z (z, 9)) is local diffeomorphism 
near any given fixed point zq G O. 

i5nc? o/ £/ie proof of Theorem U.tH Fix 6*o as above, so that 7^ —f' z (z, 9q) £ WF a (g) for all z in some 
neighborhood of f2. (Notice that g now denotes the extension by of the originally defined function 
on Q.) Let zq be a point in supp (q), where /(•, 0o)| / \ is minimal. Then —f'(z , 9 ) belongs to 

lsupp(gj ^^^^^ 

the exterior conormal cone of supp (q) at z and we get a contradiction between l|6.24l) and the fact 
that all such exterior conormal directions have to belong to WF a (q). (This is the wavefront version 
of Holmgren's uniqueness theorem, due to Hormander ([7]) and Sato-Kawai-Kashiwara (remark by 
Kashiwara in fT2]V) □ 

7 Complex geometrical optics solutions with Dirichlet data 
on part of the boundary 

In this section we prove Theorem ll.il 

We first use the Carleman estimate <|3.36ll and the Hahn-Banach theorem to construct CGO 
solutions for the conjugate operator P* = (e% Pe~%)* where * denotes the adjoint. Notice that 
P* has the same form as P v except that q is replaced by q and ip by — ip. 

Proposition 7.1 Let (p be as in fUb}) . Let v S H' 1 ^), u_ S L 2 (dQ-; {-(p' ■ v)S{dx)). Then 
3u e H°(Q) such that 



Moreover 



\\u\\ H0 + y/hWtf ■ ^)-^||an + < C^Mh-! + Vh\\(-p' ■ i>)-*v-||an_). (7.1) 

Proof: We use the Carleman estimate <|3.36[1 . Let v as in the proposition. For w £ (Hq n H 2 )(il) 
we have 

\{w\v)q + (hd„w\v-) d n_\ < \\w\\ H i\\v\\ H -i + ((-tp' ■ v)*hd v w\(-y' ■ u)~*V-) 
Therefore 

\(w,v)n+(hd u w\v-)dn_ \ < C(-i-||z;||j f -i/z||z(j|| jW i + — ^=|| ( — t^'-^) - -\/^|| ( — v 7 ' ^" ^o>^^ 1 1 )- 
Now by using H3.36[) we get 

\(w,v) a +{hd v w\v-)gaj < C(~\\v\\ B -i+-j=\\(-(f/-v)-iv- \\ga_)(\\P v w\\+Vh\\(i//-v)ihd v w\\g a+ ). 
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By the Hahn-Banach theorem, 3u 6 H°(fl), u + E L 2 (<9f2 + , ((p' ■ v) ?dS), u+ on d£l+ such that 
(w,v) a + (hd v w\v-) m _ = (P v w\u) + {hd v w\u + ) an+ , Vw e (H^ n H 2 )(n) (7.2) 

with 

\H\ H o + -j^W ■ v)~*u+\\en+ < C{-\\v\\ H -i + -^IKV ■ v)-*v-\\aa_). (7.3) 

Since P v = —h 2 A+ a first order operator, and w\ dn — we have (P v w\u) — (w\P*u) — 

h 2 (d u w\u) d n- 

Using this in (|7.2I) we obtain 

= (w\v - P*u) + h((d u w\l d n_v-) d n - (d v w\lgn + u+)an + {d v w\hu) d a) 

where lan± denotes the indicator function of dQ±. 
By varying w in (H^ n H 2 )(Q) we get 

P*u = v, hu\ m = -loci-V- + len + u+. 

which implies the proposition after replacing u_ above by —hi)—. □ 

Let W- C c?r2_ be an arbitrary strict open subset of <9f2_. We next want to modify the choice 
of u 2 in 1)5.2(1 so that u 2 \ w = 0. 

Proposition 7.2 Lei a 2 , if, tp 2 be as in Then we can construct a solution of 

Pu 2 = 0, u 2 \ w _=0 (7.4) 

of the form 

u 2 = e^ +i ^(a 2 +r 2 )+u r (7.5) 
where u r = e l ib(x; h) with b a symbol of order zero in h and 

Iml(x) = -<p{x) + k{x) (7.6) 

where k(x) ~ dist (x, 9f2_) in a neighborhood of dQ- and b has its support in that neighborhood. 
Moreover, \\r 2 \\ H o = 0(h), r 2 \ gQ _ = 0, \\((p' ■ v)~^r 2 \\ d n + = 0(hi). 

Proof: We start by constructing a WKB solution u in of 

-h 2 Au = 0,u\ dQ _ =e^+ i ^( X a 2 )\ dn _ (7.7) 

where x € C£°(dfi_), X = 1 on W-. 

We try u = e^ x 'b(x; h). The eikonal equation for I is 

((") 2 = to infinite order at dil 
= fa-tip. 

Of course g := ip 2 ~itp is a solution but we look for the second solution, corresponding to having 
u equal to a "reflected wave" . We decompose on <9f2_ 

g' = g't + sf v 
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where t denotes the tangential part and v the normal part. 
Then in order to satisfy the eikonal equation we need 

Therefore we can solve l|7.8|) to oo-order at c?Sl_ with I satisfying 

^ldn_ = 9\ an _id lJ l\ an _ — —d v g\ m _. 
By the definition of <9£1_ we have 

d„Img — —d v <~p > on 9fL_. 

Since v is the unit exterior normal we have that (|7.6|) is satisfied. 

Solving also the transport equation to oo-order, at the boundary we get a symbol b of order 
with support arbitrarily close to suppx, such that 

f -h 2 A(e%b(x;h)) = e7T0((dist ( x , dn))°° + h°°) 



= e s x«2 
an 



Our new WKB input to 1*2 will be 

(eT?a 2 - e'^b). 

Instead of l|4.18ll we get 

P(e l *a 2 -e*6) = e*fr 2 d (7.9) 

where d = 0(1) in £ 2 (f2). 

Using Proposition 17 . II we can solve 

e *Fe k (e k r 2 ) = —h d 
^\ dn _ = 

IN| H o + JhW ■ v)- i r2\\on + < j\\h 2 d\\ H -i = O(h). 

\\r 2 \\ = O(h), • i/)-i^|| 8n+ = O(V^). (7.10) 

U2 = eTZ^+i^)( a2+ r 2 )-e%b. (7.11) 
Clearly Pu 2 = 0, u 2 | an = in W— □ 
Proof of Theorem \l.l\ Let u 2 = m 2 be as in Prop 17.21 Let u\ G i? 1 (fi) solve (|5.4f) 

(A-gi)ui = 0, ui| ao = "2| an . 

By construction we have that suppit^^ n W- = 0, i = 1,2. As in Sectional let u = U\ — u%, 
q = qi — (72 ■ Then (|5.6() and (|5.7() are valid and in fact u G H 2 (il) so that the Green's formula 
l|5.8l) is also valid. Now choose v as in 15.9fl . (|5.1U|) . Then instead of (|5.11l) we get 

J a qeU^+^)( ai + n)(a 2 + r 2 )dx - J n gef ~* + i ^b(a 1 + n)dx 

= Ian + ,Jd l ,u)e-i(^(a 1 +r 1 )dS. [ " J 



with 
Thus 

Now we take 
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The second term of the LHS is what is different from l|5.11|l . Because of <|7.6[) this term goes to 
as h goes to zero, since 

\e h t+ k | = e h , 

and q, 6, ai, are bounded and || — > 0, h — > 0. Therefore we get, instead of l|5.16[l . 

Jan + , eo (^") e " * ) (^TTT) S(dx) | 2 
<^ll«i + r 1 |||n + jM<rf- 
The previous estimates imply that 

||e _ £gtt 2 ||, ||oi + ri||en + ,e = O(l). 
Consequently the RHS of 1|7.12[1 tends to as h — > and we get (|5.18|l as before, namely 

/ q{x)a 2 {x)^Jx)e lf{x) dx = 0. (7.13) 

Now the arguments of Section pimply that q = finishing the proof of Theorem ll.il □ 
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